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INTRODUCTION 

The  distribution  of  permafrost  is  controlled  by  a 
number  of  interrelated  environmental  factors  and  pro¬ 
cesses.  The  relationship  of  these  factors  and  processes 
in  Interior  Alaska,  where  permafrost  occurs  discon- 
tinuously,  is  complex  and  poorly  defined. 

A  study  was  initiated  in  1 985  by  CRREL  to  examine 
these  relationships  in  the  Caribou-Poker  Creeks  Re¬ 
search  Watershed  (Fig.  1),  situated  in  upland  terrain 
north  of  Fairbanks.  Within  this  watershed,  permafrost 
and  non-permafrost  terrain  can  be  found  having  varying 
slope,  aspect,  vegetation  and  elevation,  affording  the 
opportunity  to  analyze  the  importance  of  cover  and 
ground  conditions  in  controlling  ground  temperatures 
and  frozen  ground  distribution  (Collins  et  al.  1988). 

During  drilling  to  install  thermistor  cable  assemblies 
at  representative  sites  in  the  watershed,  an  anomalous 
unfrozen  zone  was  identified  in  a  borehole  located  on  a 
north-facing  slope.  This  unfrozen  zone  was  associated 
with  wet  conditions  that  included  the  flow  of  a  water  and 
sediment  mixture  or  slurry  to  the  ground  surface  from  a 
depth  of  about  8.5  m,  suggesting  a  possible  relationship 
to  groundwater.  Subsequent  monthly  downhole  tem¬ 
perature  measurements  (described  later)  also  show  that 
these  materials  remained  unfrozen  throughout  the  year. 

The  presence  of  unfrozen  zones  in  a  north-facing 
slope  is  considered  an  anomaly  since  the  general  con¬ 
cept  is  that  north-facing  slopes  in  Interior  Alaska  are 
usually  continuously  frozen  while  south-facing  slopes 
are  free  of  permafrost. 

Geophysical  surveys  by  impulse  radar  have  been 
used  to  resolve  the  groundwater  table  and  special  hy¬ 
drologic  features  in  areas  of  seasonally  and  perennially 
frozen  sediments  under  certain  conditions  (Annan  and 
Davis  1976,  Arcone  and  Delaney  1987,  Arcone  et  al. 
1989,  Delaney  et  al.  1990).  Under  the  hypothesis  that 
the  unfrozen  zone  was  caused  by  groundwater,  we 
undertook  an  impulse  radar  survey  of  the  site  to  deter¬ 
mine  if  the  thaw  zone  could  be  identified,  to  define  its 
distribution  and  to  determine  its  origin.  This  report 


presents  the  results  of  this  survey  in  light  of  existing 
geological  and  related  data,  and  hypothesizes  on  the 
origin  of  the  unfrozen  zone. 


BACKGROUND 
Field  site 

The  Caribou-Poker  Creeks  Research  Watershed  is 
located  48  km  north  of  Fairbanks,  Alaska,  within  the 
Yukon-Tanana  uplands  (Fig.  1).  This  region  is  char¬ 
acterized  by  rounded,  even-topped  hills  and  ridges  with 
gentle  side  slopes  leading  to  low  rel  'ef,  alluvium-floored 
valleys  (Wahrhaftig  1965).  Bedrock  in  the  watershed 
consists  mainly  of  quartz-muscovite  schist,  with  lesser 
amounts  of  garnet-mica  schist  and  micaceous  quartzite 
of  the  Fairbanks  Schist  Unit  of  the  Yukon-Tanana 
Metamorphic  Complex  (Forbes  and  Weber  1982, 
Robinson  et  al.  1990).  The  upper  1  to  20  m  of  the  rock 
is  heavily  weathered  and  covered  by  colluvium  derived 
from  it.  Loess  thinly  mantles  the  colluvium  in  some 
areas,  but  much  of  it  has  been  transported  downslope 
and  redeposited  in  valley-bottom  fans  (Rieger  et  al. 
1972). 

Vegetation  in  the  watershed  is  typical  of  an  Interior, 
upland  boreal  forest,  with  black  spruce  mainly  growing 
on  north-facing  slopes  and  poorly  drained  valley  floors, 
while  aspen,  birch,  alder  and  white  spruce  occupy  well- 
drained  soils  common  on  south-facing  slopes.  Biological 
communities  are  locally  more  complex  in  response  to 
variations  in  slope,  aspect,  drainage  and  other  factors. 

Borehole  and  radar  profile  locations 

Boreholes  were  drilled  in  April  1985  at  eight  loca¬ 
tions  on  the  southern  and  northern  faces  of  Caribou 
Peak  (elev.  773  m)  (Collins  1986).  Three  boreholes  are 
located  within  the  area  investigated  for  this  study:  K20, 
K20A  and  K2  (Fig.  2).  Borehole  K2  lies  at  the  top  of 
Caribou  Peak.  Boreholes  K20  and  K20A  lie  about  150 
m  apart  along  a  line  trending  approximately  perpen¬ 
dicular  to  the  slope  and  located  at  an  elevation  of  7 1 0  m 


Figure  1.  Location  map  showing  the  Carihou-Poker  Creeks 
Research  Watershed  within  the  Fairbanks  region  of  Interior 
Alaska. 


about  60  m  (400  m  distance)  below  K2  (Fig.  2).  The 
estimated  angle  of  slope  increases  from  2°  near  K2  to 
1 0°  at  K20-K20A  and  to  1 5°  approximately  30  m  below 
K20-K20A.  Borehole  K20  occurs  within  the  unfrozen 
zone,  while  K2  and  K20A  penetrate  into  permanently 
frozen  ground  below  the  active  layer.  Details  of  the 
borehole  logs,  site  vegetation  and  thermistor  installation 
are  presented  in  Collins  (1986)  and  Collins  et  al.  ( 1 988). 

Four  transects  for  impulse  radar  profiling  were  es¬ 
tablished  by  tape  and  compass,  with  the  lowermost 
transect  ( 1 )  located  along  a  west-east  line  between  K20 
and  K20A  that  is  perpendicular  to  the  slope  (Fig.  2). 
Transects  were  radar  profiled  on  4  April  1990,  when 
snow  ranged  in  depth  from  1.5  to  2.5  m  (Fig.  3).  Snow 
along  each  transect  was  packed  down  using  snowshoes 
prior  to  ground  profiling.  Each  transect  was  profiled 
twice  with  two  different  antenna  systems. 

Radar  instrumentation 

General  operation 

Commercially  available  short-pulse  radar  (also 
known  as  impulse,  subsurface  or  ground-penetrating 
radar)  consists  of  a  control  unit,  antennas  and  cables,  a 
magnetic  tape  recorder  and  a  power  supply.  The  control 
unit  generates  timing  signals  to  key  the  transmitter  on 


and  off  and  synchronizes  this  keying  with  the  receiver. 
It  controls  the  scan  rate  (how  fast  individual  echo  scans 
are  compiled),  the  time  range  over  which  one  wants  to 
view  the  echoes,  and  the  gain  to  be  applied  to  the  echoes. 
The  transmit  antenna  is  excited  by  a  semiconductor 
device  and  is  designed  to  radiate  a  broadband  pulse  of 
only  nanoseconds  duration.  Consequently,  the  antenna 
produces  a  low  gain  radiation  pattern.  A  separate  but 
identical  receive  antenna  is  employed  because  echoes 
can  return  from  near-surface  targets  before  the  transmit 
antenna  has  stopped  radiating.  Electronics  for  the  trans¬ 
mitter  and  receiver  are  usually  incorporated  into  the 
antenna  housing.  The  VHF-UHF  received  signals  are 
converted  by  sampling  into  an  audio  frequency  facsimile 
for  filtering,  amplification  and  recording.  Data  can  be 
simultaneously  displayed  or  quickly  played  back  in 
strip  chart  form.  Digital  recording  was  used  in  this 
study,  along  with  packaged  signal  processing  software 
that  permits  several  filtering  functions  and  data  display 
formats. 

Control  unit 

The  radar  system  used  was  manufactured  by  the 
GSSI  company  (Geophysical  Survey  Systems,  Inc., 
North  Salem,  New  Hampshire).  The  control  unit  was  a 


SIR  Model  4800  mainframe  that  triggers  pulses  at  a 
repetition  frequency  of  approximately  50  kHz  and 
compiles  the  received  pulses  into  25.6  scans/second 
(higher  or  lower  rates  are  possible).  A  variety  of  range 
gain  functions  may  be  applied  over  the  time  duration  of 
the  scans  to  suppress  the  higher  amplitude,  early  returns 
(especially  the  direct  coupling  between  transmit  and 
receive  antennas)  and  enhance  the  lower  amplitude, 
later  returns.  The  control  unit  also  controls  the  time 
range  of  the  scans,  which  can  be  varied  from  tens  to 
thousands  of  nanoseconds.  An  event  marker  is  available 
to  mark  distance  on  the  radar  record. 

Antennas 

Two  antenna  sets  (Models  3102  and  3105,  GSSI, 
Inc.)  operating  over  different  frequency  bands  were 
used  in  this  study. 

A  typical  radiated  wavelet  (Fig.  4)  for  both  antenna 
units  is  commonly  characterized  by  the  frequency 
corresponding  to  the  period  of  one  oscillation;  500  MHz 
for  the  3 1 02  radiating  in  air  and  about  300  MHz  for  the 
3105,  with  both  waveforms  having  a  bandwidth  of 
about  50%.  When  directly  coupled  to  the  packed  snow, 
these  frequencies  lowered  to  about  330  and  1 50  MHz 
respectively.  Our  use  of  two  different  frequency  bands 
was  motivated  by  our  interest  in  achieving  high  resolu¬ 


tion  near  the  surface  with  the  3 1 02,  and  deeper  penetra¬ 
tion  with  the  3105,  the  lower  frequency  of  which  would 
experience  less  attenuation  from  scattering  and  conduc¬ 
tive  absorption. 

Both  transmit  and  receive  antennas  were  housed  in 
single  units,  in  which  they  are  adjacent  and  backshielded. 
The  antennas  are  flared  dipoles  with  a  resistive  loading 
to  inhibit  current  oscillation,  and  are  excited  by  a 
transmitter  that  radiates  about  8  W  (3102)  and  5  W 
(3105)  peak  power  as  rated  by  the  manufacturer.  Sub¬ 
surface  radiation  patterns  for  either  antenna  design  are 
believed  to  be  multilobed  based  on  the  theoretical 
results  for  point  dipoles  by  Engheta  et  al.  (1982).  This 
theory  shows  a  strong  lobe  of  energy  to  be  confined 
within  the  angle  \|/  =  2sin”>(l/n)  measured  about  verti¬ 
cal,  where  n  is  the  real  part  of  the  refractive  index  of  the 
ground.  For  frozen  ground  of  high  ice  content,  n  might 
typically  be  2.3  so  that  Vj/  =  52°. 

Data  processing,  display  and  interpretation 

All  data  were  recorded  in  scans  of  512  eight-bit 
words  on  a  GSSI  DTbOOO  tape  recorder.  The  recorder 
can  be  read  onto  a  personal  computer  by  a  software 
program  (Boucher  and  Galinovsky  1989)  that  also 
processes  the  data.  The  program  is  capable  of  a  variety 
of  filter  functions,  both  horizontal  (i.e.,  over  distance  or 


Figure  4.  Typical  waveform  of  the  radar  pulse  emitted  from  either  the 
Model  3105  or  3102  antenna  units.  The  time  scale  is  arbitrary. 
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Figure  5.  Wiggle  trace  display  (left)  of  a  single  scan  and 
equivalent  graphic  record  (right)  should  these  scans 
remain  constant  with  profile  distance. 


equivalently,  many  scans)  and  vertical  (i.e.,  over  time, 
or  equivalently ,  along  one  scan).  Data  can  be  displayed 
in  either  a  wiggle  trace  (amplitude  versus  time  of  each 
scan)  or  line  scan  (intensity  as  color  or  darkness  versus 
time)  format  as  shown  in  Figure  5.  Individual  scans  are 
easily  retrievable  with  the  software.  Interpretation  of 
data  is  generally  based  on  the  simple  echo  delay  formula 

d  =  ctl2n  ( 1 ) 

where  d  is  the  depth  of  a  reflector  in  centimeters,  t  is  the 
echo  time  delay  in  nanoseconds  and  c  is  the  speed  of 
electromagnetic  waves  in  a  vacuum  (30  cm/ns).  The 
quantity  n  is  often  replaced  by  Ve  where  e  is  the  di¬ 
electric  constant.  The  factor  of  two  in  eq  1  accounts  for 
the  round  trip  propagation  path  of  the  pulse.  Equation  I 
applies  only  to  reflections  emanating  from  horizontally 
flat  interfaces  of  several  in-situ  wavelengths  extent,  or 
to  scattering  from  point  sources.  It  can  be  applied  to 
several  layers  successively  if  n  is  known  for  each  layer 
and  the  time  delays  to  each  layer  interface  are  easily 
picked  off  the  record. 

The  dielectric  constant  of  frozen  or  unfrozen  earth 
materials  as  a  function  of  water  content  in  the  100-  to 
10(X)-MHz  range  has  not  been  extensively  investigated 
in  comparison  to  the  wide  variety  of  material  types  that 
exist.  Generally,  competent  crystalline  rocks  can  range 
from  E  =  4  (quartzite)  to  6  (felsic  granites)  to  9  (marble 
or  basalts;  e.g.,  Campbell  and  Ulrichs  1969)  with  little 
dependence  on  water  content  as  it  is  generally  less  than 
\%  of  the  volume  of  the  material.  There  are  no  data 
available  for  deeply  weathered  bedrock,  but  it  should  be 


similar  to  the  properties  of  soil  for  which  there  are 
extensive  data,  both  laboratory  and  field.  These  inves¬ 
tigations  have  shown  that  e  of  unfrozen  wet  soils  (sands, 
silts,  clays,  loams)  will  monotonically  increase  with 
increasing  water  content,  and  will  range  from  about  2.6 
when  dry,  to  about  36  at  saturation  (e.g.,  Topp  et  al. 
1980),  with  dielectric  loss  increasing  in  similar  pro¬ 
portion.  In  contrast,  the  e  of  frozen  soils  will  increase 
only  until  about  40%  water  content  is  reached,  beyond 
which  ice  begins  to  replace  unfrozen  adsorbed  water 
(which  is  significant  in  silts  and  clays)  and  thus  decreases 
e  (e.g.,  Arcone  and  Delaney  1989,  Delaney  and  Arcone 
1984).  Typically,  e  for  the  high  ice  content  sands,  silts 
and  gravels  of  Interior  Alaska  ranges  between  4  and  6. 
Massive  ice,  of  course,  has  an  e  =  3.2-3.5  and  a  dry 
snowpack  about  1.2-1. 4  (Gumming  1952)  for  typical 
densities  of  a  seasonal  ground  cover. 

Reflections  that  strongly  slope  or  curve  with  the 
horizontal  profiling  distance  must  be  “migrated,”  a 
complicated  process  of  diffraction  summation  that  de¬ 
rives  the  correct  orientation  of  reflectors  from  the  radar 
time-distance  data.  Migration  is  a  straightforward  pro¬ 
cedure  when  there  is  only  one  layer  of  interest.  For  a 
more  complicated  layer  structure  in  which  n  may  vary 
continuously  within  the  space  of  a  few  in-situ  wave¬ 
lengths,  exact  migration  is  practically  hopeless. 
Therefore,  migration  of  sloping  reflectors  in  this  study 
will  assume  a  simplified  structure  for  n  based  upon  data 
available  in  the  above  cited  references. 

Borehole  geology 

The  stratigraphic  log  of  borehole  K20  is  typical  of 
the  watershed  uplands  (Fig.  6).  A  0.15-  to  0.1 8-m 
organic  mat  overlies  a  relatively  thick  (3.0  m)  layer  of 
gravelly  colluvium  containing  very  angular  rock  frag¬ 
ments,  which  in  turn  overlies  a  0.4-m  layer  of  gravelly 
silt.  This  horizon  overlies  quartz-muscovite  schist 
bedrock,  which  is  weathered  to  a  depth  of  about  1 3.9  m 
(Collins  1986).  Moderately  weathered  bedrock  grades 
progressively  below  a  depth  of  about  10  m  to 
unweathered  bedrock  at  13.9  m  where  auger  refusal 
occurred.  It  should  be  noted  that  because  of  equipment 
limitations,  the  drill  was  unable  to  penetrate  into  sound 
bedrock.  This  limitation,  however,  clearly  defines  the 
depth  of  weathering  at  each  site. 

The  ground  was  frozen  to  3.0  m  depth,  where  an 
unfrozen  layer  of  gravelly  silt  was  encountered.  This 
localized  layer  remained  unfrozen  until  May  1989, 
when  measurements  indicated  that  ground  temperatures 
had  gradually  decreased  over  several  months  to  below 
0°C.  The  gravelly  silt  horizon  was  therefore  classified 
as  frozen  during  this  radar  study,  but  not  so  strongly  as 
to  seriously  depress  the  dielectric  constant  of  the  material, 
as  will  be  discus.sed  later.  From  3.4  to  about  6.0  m,  the 
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Figure  6.  Summary  of  geologic  logs  for  boreholes  K20.  K20A  ami  K-2. 


weathered  bedrock  wa.s  dr>'  and  frozen.  Cuttings  from 
depths  of  6.0  to  8.5  m  appeared  moist,  but  it  was  un¬ 
certain  if  the  bedrock  was  frozen  or  unfrozen.  Subse¬ 
quent  temperature  measurements  suggest  that  the  bed¬ 
rock  at  these  depths  is  frozen  at  6.0  m,  but  warming  to 
slightly  above  0°C  by  the  8.5-m  depth.  Extremely  wet, 
unfrozen  conditions  occurred  below  8.5  m.  Water  con¬ 
taining  cuttings  and  weathered  materials  flowed  to  the 
ground  surface  as  a  slurry  from  a  depth  of  about  8.5  to 
9.0  m  during  the  drilling.  After  1 2  hours,  water  associated 
with  this  zone  had  risen  to  5  m  depth  in  the  borehole. 

Boreholes  K20A  and  K2  were  frozen  through  their 
approximately  5-m  lengths  in  1985  (.Collins  1986). 
Unweathered  quartz-muscovite  senist  occurred  at  a 
shallower  depth  in  K20A  (5.3  m)  than  in  K20.  and  the 
overlying  colluvium  appeared  to  grade  into  weathered 
bedrock  from  about  2. 1  to  about  3.5  m  (Fig.  6).  Borehole 
K2  consisted  of  a  much  thinner  (-1.5  m)  frozen  collu¬ 
vial  layer  overlying  slightly  weathered  quartz-musco¬ 
vite  schist;  auger  refusal  occurred  at  a  depth  of  4.6  m. 


RESULTS  AND  DISCUSSION 
Radar  profiles 

Seven  radar  profiles  at  presented  in  wiggle  trace 
format,  three  for  the  330-MHz  3102  antenna  unit  and 
four  for  the  150-MHz  3105  unit.  The  extra  3105  profile 
was  run  perpendicular  to  transects  1  and  3.  On  each 
profile,  the  horizontal  direction  is  distance  along  the 
transects  and  the  vertical  direction  is  echo  time  of 
return.  Each  profile  was  processed  with  a  horizontal 
high-pa.ss,  “background  removal”  filter  to  eliminate 
noise  bands,  and  a  vertical  high-pass  filter  to  eliminate 
variations  in  overall  signal  level  from  scan  to  scan 
caused  by  changes  in  coupling  between  antenna  and 
snow  surface  because  of  uneven  towing  of  the  antenna. 
All  profiles  were  also  horizontally  justified  in  that  the 
number  of  scans  between  distance  marks  was  equalized. 

The  events  on  the  records  are  classified  either  as 
horizons,  which  are  generally  horizontal  and  extend 
across  the  entire  record,  or  as  dipping  reflectors  that 
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might  either  be  real,  or  apparent,  as  would  be  caused  by 
air  or  snow  refractions.  The  horizons  are  composed  of 
intermittent  reflections  mat  occur  frequently  enough  to 
indicate  a  vertical  transition  in  material  properties, 
whereas  the  reflections  from  the  dipping  structures  are 
much  stronger  and  more  consistent.  The  150-MHz 
profiles  have  longer  time  scales  (as  allowed  by  the  data 
digitization  rate)  to  look  for  deeper  reflections.  There¬ 
fore,  similar  events  in  both  profiles  only  appear  to  occur 
at  different  depths  or  have  different  time-distance  slopes. 

Figure  7  shows  the  1 50-  and  330-MHz  profiles  along 
the  190-m  transect  1 ,  which  passed  directly  adjacent  to 
boreholes  K20  and  K20A  (Fig.  2).  The  330-MHz  record 
has  a  shorter  time  duration  and  therefore  the  events  are 
better  separated  than  in  the  ISO-MHz  record.  The  first 
horizon  (number  1)  is  the  direct  coupling  between  the 
transmit  and  receive  antennas  and  serves  as  a  0-ns 
reference  for  timing  the  later  events.  Horizon  2  is  the 
base  of  the  snow  cover  and  generally  occurs  at  about  5- 
to  7-ns  delay.  If  the  trampled,  seasoned  snowpack  has  a 
density  of  about  0.35,  the  corresponding  dielectric 
constant  is  about  1.4  (Camming  1952),  which  translates 
the  time  delays  to  snow  depths  of  89-124  cm. 

Horizon  3  is  far  less  consistent  and  appears  mainly  as 
a  darkening  of  the  record;  it  is  more  obvious  in  the  330- 
MHz  profiles  of  transects  2  and  3  shown  later.  We 
interpret  this  horizon  as  a  transition  within  the  colluvium 
and  not  between  the  colluvium  and  weathered  bedrock, 
which  is  far  too  deep  (3.0-3.4  m)  to  give  a  meaningful 
dielectric  constant  for  the  time  delay  of  about  25  ns  near 
K20.  Using  a  dielectric  constant  of  about  4.8,  as  mea¬ 
sured  in  gravelly  alluvium  (Arcone  and  Delaney  1989) 
of  similar  water  content  (5-30%  by  weight)  and  tem¬ 
perature  (0  to  -1°C),  places  this  horizon  at  a  depth  of 
about  1.7  m  below  the  ground  surface.  Therefore,  this 
horizon  more  likely  indicates  a  seasonally  derived 
transition  from  low  (about  3%  or  less)  to  high  (about  3- 
6%)  volumetric  unfrozen  water  content,  as  is  known  to 
occur  in  fine-grained  sediments  (Hoekstra  et  al.  1975); 
a  few  percent  change  in  unfrozen  water  content  has  been 
shown  both  theoretically  and  experimentally  by  Delaney 
and  Arcone  ( 1 984)  to  approximately  double  the  dielectric 
constant.  This  horizon  is  deeper  further  across  the 
record,  reaching  a  depth  of  about  2. 1  m  near  K20A.  The 
transition  in  unfrozen  water  content  cannot  be  deduced 
from  the  drill  log  that  documents  only  total  water 
content  in  cuttings  at  specific  locations.  There  is  no 
distinct  horizon  forthecolluvium-to- weathered-bedrock 
transition,  which  suggests  that  the  weathered  bedrock 
has  a  dielectric  constant  similar  to  that  of  the  marginally 
frozen,  lower  layer  of  colluvium. 

The  most  conspicuous  features  in  both  the  profiles  of 
Figure  7  are  the  dipping  reflectors.  Those  that  dip  most 
steeply  eithe;  to  the  east  or  west  (e.g.,  event  4)  are 


refracted  air  waves  traveling  along  the  air/snow  bound¬ 
ary  and  reflected  back  to  the  antenna.  This  source  is  de¬ 
duced  from  the  time/distance  slope,  which  gives  a  one¬ 
way  wave  velocity  very  close  to  30  cm/ns,  the  speed  of 
radio  waves  in  a  vacuum.  The  stronger  reflections  (e.g., 
events  5  and  6),  which  all  dip  to  the  east,  are  believed  to 
be  legitimate  linear  features  that  are  probably  water 
within  fractures  in  the  bedrock,  which  drill  logs  suggest 
is  moderately  to  slightly  weathered.  Three  primary 
planes  appear  evenly  spaced  at  about  55  to  60  m.  The 
steepest  planar  feature  at  0  to  8  m  along  profile  1  is  a 
multiple  of  the  feature  to  its  right.  In  either  profile,  the 
strongest  part  of  these  events  occurs  at  a  time  delay  of 
about  120  ns  between  stations  20  and  30  m  (event  6).  At 
this  point  there  is  also  the  peak  of  a  hyperbola,  which 
indicates  a  strong,  localized  concentration  of  water,  as 
might  occur  in  an  open  fracture. 

Since  fractures  would  be  characteristic  of  only  the 
bedrock,  we  may  assume  that  the  start  of  these  fea¬ 
tures  indicates  the  colluvium/weathered  bedrock  inter¬ 
face.  These  features  seem  to  originate  at  progressively 
greaterdepths  from  the  surface  when  moving  from  west 
to  east  across  Figure  7,  but  the  exact  depths  cannot  be 
determined  since  we  have  no  independent  measure  of 
the  dielectric  constant  of  the  entire  column  of  collu¬ 
vium.  Nor  is  the  dielectric  constant  of  the  weathered 
bedrock  known,  which  is  needed  to  determine  the 
dip  angle  of  the  fractures.  However,  this  latter  con¬ 
stant  need  only  be  roughly  estimated  to  find  an  approxi¬ 
mate  dip  angle. 

In  a  homogeneous  medium,  the  apparent  dip  angle  a 
in  nanoseconds  per  meter  measured  from  horizontal  in 
a  time-distance  plot  is  related  to  the  true  dip  angle  P 
in  a  depth-distance  plot  through  the  migration 
relationship 

P  =  sin~'  (ctana/2«).  (2) 

This  angle  will  not  change  on  the  radar  record  if  there  is 
an  overburden  of  constant  depth  and  dielectric  profile 
situated  above  the  medium  containing  the  fracture.  If 
there  were  a  direct  correspondence  in  scale  of  1  m  =  1 
ns,  the  apparent  fracture  angle  a  in  Figure  7  would  be 
approximately  64°.  As  there  is  no  evidence  of  a  strong 
bedrock-colluvium  reflection,  and  the  reflection  within 
the  colluvium  is  weak  at  best,  n  for  the  weathered 
bedrock  probably  varies  only  between  about  2  and  3, 
confining  P  to  between  6  and  9°. 

Borehole  K20  encountered  water  at  a  depth  of  about 
8.5  to  9.0  m.  Using  n  =  2  for  the  weathered  bedrock, 
which  gives  a  true  dip  angle  p  =  9°,  would  place  the  radar 
planar  feature  at  about  8-9  m  below  the  ground  surface 
at  K20  and,  thereby,  links  the  wet,  unfrozen  zone 
directly  to  these  dipping  reflectors.  The  similarity  in 
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Figure  7.  Processed  radar  profiles  of  transect  1 .  Boreholes  K20  and  K20A  are  located.  Processing  and  explanations  of  the  numbered  events  are  described  in  the  text . 
The  strong,  relatively  planar  signatures  dipping  to  the  right  (east)  are  interpreted  as  water-fdledfractures  within  the  bedrock,  a.  330-MHz  3102  antenna,  b.  150-MHz 
3105  antenna. 
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Figure  8. 330-MHz  profile  of  transect  I ,  with  additional  stacking 
and  fdtering  (left),  to  better  illustrate  features  of  interest  as 
interpreted  on  the  right. 


returns  of  the  other  planes  indicates  that  they  are  also 
probably  filled  with  water  at  depth.  Downhole  tempera¬ 
ture  data  provide  further  evidence  on  the  nature  of  this 
zone  and  are  discussed  later. 

Figure  8a  is  a  highly  processed  version  of  Figure  7a, 
and  its  interpretation,  which  summarizes  many  of  the 
ideas  presented  above,  is  given  in  Figure  8b.  The  data  of 
Figure  8a  were  first  stacked  six-fold  and  then  processed 
using  a  Hilbert  magnitude  transform,  which  is,  basi¬ 
cally,  an  energy  envelope  detection  scheme.  The  pur¬ 
pose  of  the  stacking  is  to  compress  the  data  horizontally 
to  bring  out  the  horizons.  Figure  8b  interprets  the 
horizons  of  Figure  8a  as  the  base  of  the  snow  layer,  a 
transition  from  low  to  high  unfrozen  water  content 
within  the  colluvium  (horizon  3,  Fig.  7a),  and  the 
weathered  bedrock  as  interpreted  from  the  start  of  the 
fracture  events.  The  low-to-high  unfrozen  water  tran¬ 
sition  is  seen  to  rise  in  the  middle  of  the  record  and 
deepen  towards  the  east.  This  rise  in  the  center  of  the 
record  implies  that  borehole  K20  is  not  located  where 
the  warmest  zone  comes  closest  to  the  surface  along  this 
profile. 


Two  additional  transects  (transects  3,  4;  Fig.  2) 
were  profiled  to  determine  if  the  planar  features  at  tran¬ 
sect  1  were  highly  localized  or  more  continuous.  Each 
of  the  two  profiles  that  lie  parallel  to  transect  1  exhibit 
several  planar  features,  with  the  strongest  continuous 
reflection  correlating  precisely  with  the  structure  in¬ 
tercepted  by  borehole  K20  (Fig.  9).  Profile  3,  located  60 
m  upslope  of  profile  1 ,  also  has  a  strong  return  at  depth 
(about  1 10-130  ns)  that  appears  to  correlate  with  the 
middle  signature  (60-80  m)  on  profile  1.  The  planar 
signature  furthest  east  on  profile  1  (about  122  m)  is  not 
as  strong  a  reflection  and  also  does  not  appear  to  have 
been  encountered  on  transect  3.  This  absence  may 
simply  have  resulted  from  transect  3  being  of  insuffi¬ 
cient  length  to  have  passed  over  it.  Similarly,  on  transect 
4,  the  profile  reveals  a  single,  strong  planar  reflector  that 
we  have  correlated  with  the  middle  signature  on  profile  1 . 
The  partial  feature  on  the  north  end  of  transect  2,  which 
is  perpendicular  to  transect  1,  may  correlate  with  a 
fourth,  only  partial  return  at  approximately  105  m  dis¬ 
tance  and  1 10  ns  depth  on  line  1 ,  but  this  is  not  certain; 
no  correlative  signature  occurs  on  either  profiles  3  or  4. 
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Fi'^ure  10.  Spatial  distribution  of  fractures  relative  to  radar 
transects  and  boreholes.  Fractures  trend  generally  N20°E  with  an 
easterly  dip  of  about  9°. 


The  spatial  continuity  and  fonn  of  the  dipping  reflec¬ 
tors  are  evidence  of  these  features  being  fractures,  or 
possibly  joints,  owing  to  their  relatively  consistent 
spacing,  within  the  bedrock.  Forbes  and  Weber  (1982) 
indicate  that  bedrock  in  the  crest  area  forms  the  east- 
west  trending  axis  of  an  overturned  or  recumbant  fold. 
The  fractures  apparently  occuron  its  northern  overturned 
limb  and  may  be  associated  with  this  deformation. 

The  distribution  of  the  fractures  and  their  apparent 
correlation  from  profile  to  profile  at  a  given  depth  (in  ns) 
indicate  that  they  trend  roughly  N20°E,  upslope  and 
downslope,  at  this  site.  The  apparent  convergence  of  the 
two  best  defined  fractures  may  be  real  or  be  ascribable 
to  ground  and  bedrock  surface  relief  for  which  the 
trends  are  not  corrected  (Fig.  10).  It  is  also  possible  that 
fractures  within  the  bedrock  of  the  surveyed  area  are 
more  numerous  than  displayed  by  the  radar.  This  pos¬ 
sibility  is  suggested  by  the  apparent  lack  of  continuity 
of  certain  signatures  from  profile  to  profile.  Unless  a 
change  in  the  dielectric  constant  of  the  materials  results 
from  the  presence  of  water  in  the  fractures,  the  fractures 
would  probably  be  indistinguishable  from  the  sur¬ 
rounding  weathered  bedrock. 

The  continuity  of  the  fractures  upslope  and  downslope 


suggests  that  groundwater  is  flowing  from  a  source  that 
is  higher  in  elevation  than  the  area  surveyed  in  this 
study,  or  alternatively,  may  have  an  aquifer  source  at 
depth  in  the  bedrock .  Aerial  photographs  of  Caribou  Peak 
during  summer  (see  Fig.  1 2  later)  show  both  outcrop 
and  vegetation  patterns  with  an  approximately  parallel 
north-south  trend,  suggesting  that  a  relationship  exists 
between  drainage  from  the  fractures  and  the  vegetation 
pattern;  this  relationship  will  be  considered  later. 

Groundwater  in  unfrozen  schistose  bedrock  else¬ 
where  in  this  region  commonly  occurs  within  structural 
features.  Cederstrom  ( 1 963)  and  P6we  and  Bell  ( 1 975), 
for  example,  concluded  that  its  occurrence  in  the 
Fairbanks  area  uplands  was  commonly  controlled  by 
fractures.  Joints  and  fractured  quartz  veins,  in  some 
cases  being  perched  above  permafrost.  Water-bearing 
structures  of  this  type  are  generally  random  in  location 
and  frequency  of  occurrence,  and  therefore  difficult  to 
predict.  The  radar  method  now  provides  one  technique 
for  better  defining  this  random  distribution. 

Downhole  temperatures 

Ground  temperatures  have  been  measured  continu¬ 
ously  in  each  borehole  since  1986  (Collins  et  al.  1988). 
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For  these  measurements  the  resistances  of  individually 
calibrated,  precision  thermistors  are  read  directly  and 
converted  to  temperature  using  the  modified  Steinhart 
equation.  Measurements  are  accurate  to  ±0.05°C,  with 
a  system  uncertainty  estimated  at  ±0.03°C  and  a  preci¬ 
sion  of  ±0.0 1°C.  Data  from  K20A,  K20  and  K2  for  the 
period  of  January  1987  through  December  1989  are 


Figure  1 1 .  Downhole  temperatures  for  1989.  Seasonal 
variations  are  attenuated  in  K20,  exhibiting  uniform, 
relatively  stable  temperatures,  with  the  unfrozen  hori¬ 
zon  occurring  at  about  8.0  to  13.5  m.  K20A  and  K2  show 
more  typical  seasonal  responses  within  about  5  m  of  the 
surface. 


summarized  i  n  Table  1 .  S  ummary  graphs  of  the  downhole 
temperatures  for  1989  are  shown  in  Figure  1 1. 

Although  the  ground  temperature  data  are  limited  in 
depth  and  quantity,  qualitative  assessments  of  their 
bearing  on  the  origin  of  the  unfrozen  zone  can  be  made. 
Three  aspects  of  the  temperature  data  from  K20  apptear 
anomalous  when  compared  with  the  data  from  nearby 
boreholes  K20A  and  K2  as  well  as  from  other  boreholes 
in  the  watershed  presented  in  Collins  et  al.  (1988). 

The  three  anomalous  aspects  of  the  temperature  data 
are  the  range  and  magnitude,  the  gradient,  and  seasonal 
variations.  The  data  suggest  that  an  influx  of  heat  not 
correlated  with  seasonal  surface  warmings  or  with  the 
geothermal  heat  flux  is  affecting  the  ground  tempera¬ 
tures  at  K20. 

Seasonal  variations  in  ground  temperatures  are  highly 
attenuated  within  K20.  Normally,  seasonal  variations  in 
air  and  surface  temperatures  cause  subsurface  fluctua¬ 
tions  in  temperature  to  a  depth  of  20  m,  where  changes 
are  about  0.0 1  °C  with  a  time  lag  of  1  year  (e.g..  Gold  and 
Lachenbruch  1973).  More  rapid  responses  of  larger 
magnitude  occur  within  the  upper  10  m. 

Commonly,  even  though  there  may  be  small  tem¬ 
perature  variations  to  20  m  depth,  the  depth  of  zero 
annual  amplitude  can  be  assumed  to  be  at  about  10 
m,  where  the  annual  variation  is  less  than  0.1  °C,  which 
is  less  than  1%  of  the  surface  temperature  annual 
amplitude  (Lunardini  1981).  The  ground  temperature 

at  depth  will  generally  be  less  than  the  mean 
annual  ground  surface  temperature  Tg.  Thermal  dis¬ 
turbances  will  modify  this  relationship  (Brown  and 
Pewd  1973).  At  the  crest  of  Caribou  Peak,  climatic  data 
indicate  that  the  Mean  Annual  AirTemperature  (MAAT) 
is  about  -2.9°C  with  a  Mean  Annual  Ground  Tempera¬ 
ture  (MAGT)  of- 1.0°C  (Haugen  et  al.  1982,  Collins  et 
al.  1988).  At  K20A,  the  mean  annual  ground  tempera¬ 
ture  is  estimated  at  -1 .3'’C  (Collins  et  al.  1988).  Bore¬ 
hole  data  suggest  the  temperature  at  10  m  depth  is  about 
- 1 . 1  °C.  The  MAAT  for  the  K20A  site  may  be  estimated 
using  the  empirical  relationship  between  MAAT  and 
MAGT  of  Haugen  et  al.  (1983) 

MAGT  =  0.7  MAAT  +  1 .25  (3) 

where  MAGT  is  approximately  the  same  as  T^,  based 
upon  a  number  of  borehole  measurements  throughout 
Interior  Alaska.  Using  a  Tj  of -1.1®C  for  MAGT,  we 
estimate  the  MAAT  value  to  be  -3.6°C. 

Given  the  proximity  of  holes  K20A  and  K20,  and  the 
sites’  general  similarity  in  elevation,  aspect,  slope  angle 
and  vegetation,  the  MAAT  and  ground  temperatures 
and  their  variations  should  be  similar  in  both  holes. 
Measurements  in  hole  K20,  however,  indicate  relatively 
high  ground  temperatures  and  seasonal  variations  that 
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Table  1.  Ground  temperature  data  from  thermistor  strings  installed  in  boreholes  K20,  K20A  and  K2. 
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♦Cable  disturbed  by  black  bear.  Temperature  data  from  July  through  Nov.  reflect  disturbance. 


Table  1  (Cont’d).  Ground  temperature  data  from  thermistor  strings  installed  in  boreholes  K20,  K20A,  and  K2. 
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are  highly  attenuated  within  3.0  m  of  the  surface  (Fig. 
1 1 ).  Below  this  depth,  ground  temperatures  do  not  show 
a  seasonal  response.  Values  increase  to  slightly  above 
0°C  with  little  systematic  variation  from  month  to 
month.  Thermistor  measurements  at  the  9.0-,  1 2.0-  and 
13.4-m  depths  indicate  an  average  annual  value  of 
0.05°C,  significantly  higher  than  anticipated  from  the 
estimated  MAAT  and  MAGT  values  for  this  location. 
Temperatures  at  this  depth  clearly  do  not  simply  reflect 
the  effects  of  local  climatic  conditions  on  the  regional 
geothermal  heat  flux. 

Temperatures  within  permafrost  can  normally  be 
defined  using  steady  state,  one-dimensional  heat  con¬ 
duction  theory.  According  to  the  Fourier  heat  conduc¬ 
tion  law  (e.g..  Gold  and  Lachenbuch  1973) 

q  =  (4) 

where  Gp  is  the  thermal  gradient  of  layer  n,  is  the 
thermal  conductivity  of  layer  /;.  and  q  is  the  rate  of  heat 
conducted  per  unit  area.  The  thermal  gradient  is  simply 
defined  as 


where  t  is  temperature  and  r  is  depth. 

Locally,  the  value  of  q  can  be  assumed  to  be  identical 
in  each  borehole.  Because  subsurface  materials  at  K20 
and  K20A  are  basically  identical  (although  the  thick¬ 
ness  of  the  colluvium  is  less  at  K20A),  the  thermal 
conductivity  of  the  stratigraphic  sequence  can  be  as¬ 
sumed  to  be  similar  as  well.  Therefore,  we  also  assume 
that  the  temperature  gradient  and  range  in  subsurface 
temperatures  at  K20  should  be  the  same  overall  as  in 
K20A  unless  modified  by  other  external  factors. 

Comparison  of  the  K20A  and  K20  data  shows  no 
similarity  to  the  temperatures  to  the  5-m  depth  nor 
projected  temperatures  for  K20A  from  the  5-m  depth  to 
the  equivalent  13-m  depth  of  K20.  Projected  ground 
surface  temperatures  are  also  dissimilar.  A  mean  annual 
surface  temperature  projected  from  the  borehole  data  of 
K20  would  be  much  higher  than  expected,  ranging  from 
-0.15  to  -0.2°C  for  1987  to  1989.  These  data  also 
indicate  a  modified  heat  flux  at  depth. 

The  relatively  constant  temperatures  and  their  lim¬ 
ited  range  suggest  that  heat  conducted  by  groundwater 
is  locally  disturbing  the  thermal  regime  in  K20  and 
maintaining  temperatures  above  the  freezing  point.  In 
the  3  years  of  record,  the  uniform  temperatures  slightly 
increased  (0.04  to  0. 1 4°C  range)  in  the  summer  months 
of  May  through  September  of  1988  and  1989,  particu¬ 
larly  in  the  months  of  June  and  July,  in  the  unfrozen 
zone  (T able  1 ).  Ground  temperatures  above  this  zone  do 
not  show  a  corresponding  increase,  and,  in  particular 


from  about  the  3. 5-m  depth  to  the  9.0-m  depth,  show 
little  variation.  In  addition  the  temperature  increases 
take  place  too  rapidly  to  be  a  response  to  seasonal 
surface  warming.  In  the  remaining  months  of  1988  and 
1989.  as  well  as  all  of  1987,  there  were  virtually  no 
variations  in  ground  temperatures  (<0.03°C)  below  3  m 
depth.  Records  from  the  nearest  climatic  station  at 
Fairbanks  indicate  that  May  through  September  are 
typically  the  wettest  and  warmest  months  in  this  region 
(NOAA  1989).  Both  1988  and  1989  had  comparable 
average  air  temperatures  during  the  summer  months, 
but  normal  or  above  average  precipitation  in  May  and, 
in  particular,  June  when  2.26  and  2.53  in.  (57.4  and  64.3 
mm)  of  precipitation,  respectively,  fell  (0.87  and  1.14 
in.  [22.1  and  28.0  mm]  above  normal).  During  1987,  air 
temperatures  were  comparable,  but  precipitation  was 
below  normal  each  month,  with  only  1 .02  in.  (25.9  mm) 
in  June  and  6.43  in.  ( 1 63.3  mm)  for  the  year,  about  5  in. 
(127  mm)  below  normal. 

Thus,  there  is  an  apparent  relationship  between  the 
amount  of  precipitation  and  the  slight  increase  in  ground 
temperatures  during  summer,  suggesting  that  the  increase 
may  have  resulted  from  infiltration  or  recharge  of  the 
groundwater  in  the  fracture  during  May  and  June. 
Relatively  constant  or  slightly  decreasing  temperatures 
would  be  expected  during  the  winter  season  when 
seasonally  frozen  ground  and  low  air  temperatures 
would  severely  restrict  or  eliminate  groundwater  re¬ 
charge.  If  groundwater  recharge  was  low  prior  to 
freezeup,  a  progressive  drop  in  groundwater  discharge 
in  the  winter  months  could  result  in  a  reduction  in  heat 
flux  and  a  decrease  in  the  temperature  in  the  outer  fringe 
of  the  thawed  zone  in  the  fall  and  winter.  Thus,  although 
the  range  in  variations  in  the  ground  temperatures  fall 
within  the  accuracy  of  the  thermistors,  the  consistent 
increase  in  the  absolute  temperatures  and  apparent 
association  with  summer  precipitation  lend  support  to 
the  idea  of  a  groundwater  source  for  the  continuously 
unfrozen  conditions. 

To  provide  some  verification  of  the  near  0°C  tem¬ 
perature  measurements,  samples  of  the  water-sediment 
mixture  from  borehole  K20  and  of  leachate  extracted 
from  the  same  sample  were  analyzed  to  determine 
whether  its  freezing  point  was  0°C  or  depre.-.sed  by  the 
presence  of  salts  to  a  more  negative  value.  This  proce¬ 
dure  was  considered  important  in  lieu  of  the  accuracy  of 
the  temperature  measurements  to  determine  if  the  water 
could  be  at  subfreezing  temperatures  and  thus  static 
within  the  permafrost. 

A  sample  from  the  water-sediment  mixture  was 
placed  in  a  controlled  temperature  cell  and  gradually 
cooled  at  0.0 1  "C  increments.  Unfrozen  water  content  was 
monitored  using  the  NMR  technique  (Tice  et  al.  1978).  It 
began  freezing  at  0°C  with  pore  ice  formation  initiated. 
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Figure  12.  Aerial  photograph  of  the  Caribou  Peak  area  on  12  July  1 977,  showing  the 
approximately  parallel  linear  arrangement  to  the  vegetation.  The  darker  lineations 
result  from  an  increased  cover  composed  of  alder,  willow  and  taller  black  spruce. 
The  fractures  identified  by  radar  profiling  (dashed  lines)  appear  to  intersect  the 
interface  between  the  weathered  rock  and  the  colluvium. 


A  leachate  sample  was  also  analyzed  by  ion  chroma¬ 
tography  to  examine  the  presence  of  principal  salt- 
forming  cations.  This  analysis  indicated  quantities  esti¬ 
mated  at  less  than  several  parts  per  million,  with  the 
probability  that  there  was  virtually  no  freezing  point 
depression. 

These  analyses  support  the  credibility  of  the  posi¬ 
tive,  near  0°C  temperature  measurements,  and  that  the 
groundwater  must  be  continuously  replenished  from  a 
fresh  water  source. 

Thus,  the  generally  stable  and  positive  near  0°C 
temperatures  within  the  fracture  zone,  the  highly  at¬ 
tenuated  seasonal  temperature  variations,  the  apparent 
lack  of  correlation  in  ground  temperatures  to  local 
climatic  conditions,  and  the  steep  temperature  gradient 
are  evidence  of  a  thermal  disturbance  caused  by  heat 
conducted  by  groundwater  flowing  within  the  fractures 
throughout  the  year. 


Vegetation  patterns  and  fracture  locations 

Changes  in  the  vegetation  covering  the  northern 
slope  of  Caribou  Peak  form  a  mostly  1  inear  pattern  (Fig. 
1 2).  Vegetative  lineations,  which  trend  generally  north- 
south  and  lie  approximately  parallel  to  one  another, 
reflect  shifts  in  the  composition  and  percentage  of 
cover.  While  the  majority  of  the  northern  face  is  covered 
by  an  open  black  spruce  forest  with  scattered  alder,  the 
vegetation  in  the  lineations  is  a  denser  cover  of  taller 
black  spruce,  alder  and  willow  These  differences  in  the 
vegetation  indicate  slightly  warmer  and  wetter  condi¬ 
tions  here  than  in  the  areas  immediately  adjacent  to 
them. 

In  addition,  the  vegetative  pattern  can  be  differenti¬ 
ated  into  those  lineations  that  are  discontinuous  and 
generally  disappear  within  the  upper  two-thirds  of  the 
slope  and  those  that  are  continuous  downslope  and 
located  within  gullies  incised  into  the  bedrock  of  about 


18 


the  lower  one-third  of  the  northern  face  (Fig.  1 2).  These 
gullies  commonly  converge,  forming  seasonal  tributar¬ 
ies  to  Poker  Creek.  Within  the  upper  two-thirds  of  the 
slope,  neither  of  these  groups  of  lineations  are  associ¬ 
ated  with  such  changes  in  surface  relief. 

The  cause  of  the  vegetative  pattern  appears  to  be 
related  to  the  fractures  within  the  bedrock.  A  plot  of  the 
approximate  positions  of  radar  profiles  on  an  aerial 
photograph  (Fig.  12)  suggests  that  the  water-bearing 
fractures  intersect  the  interface  between  the  bedrock 
and  colluvial  overburden  beneath  or  adjacent  to  the 
vegetative  lineations  (Fig.  1 3).  Thus,  the  thicker  vegeta¬ 
tion  apparently  occurs  where  higher  ground  tempera¬ 
tures  and  increased  availability  of  moisture  result  from 
groundwater  discharge  through  the  fractures. 

The  association  of  the  continuous  lineations  with  the 
intersecting  gullies  suggests  that  a  seasonal  drainage 
(there  are  no  perennial  streams  here)  results  from  sea¬ 
sonal  increases  in  groundwater  discharge  from  the 
fractures,  coupled  with  surface  runoff  from  snowmelt 
or  major  rainfall  events,  or  both.  Deeper  weathering  of 
bedrock  in  the  fracture  zones  may  be  in  part  responsible 
for  gully  development  within  the  lower  slope,  as  well  as 
for  concentrating  seasonal  surface  runoff  in  the  upper 
reaches  of  these  lineations.  The  discontinuous  lineations 
imply  enhanced  temperatures  and  moisture  conditions 
attributable  to  the  presence  of  fractures  or  joints,  but 
limited  surface  runoff  or  an  absence  of  groundwater 
flow.  The  occurrence  of  other  continuous  vegetative 
lineations  elsewhere  on  this  northern  slope  of  Caribou 
Peak  as  well  as  others  in  the  watershed  implies  that 
additional  water-bearing  fractures  and  localized  thaw 
zones  may  exist  in  this  region. 

Groundwater  origin 

The  source  of  water  flowing  within  the  fractures  is 
somewhat  problematic,  but  limited  by  the  local  geology 
and  site  geography.  Both  the  trend  of  the  fractures 
defined  by  the  radar  profiles  and  their  apparent  associa¬ 


tion  with  the  vegetative  patterns  across  the  slope  indi¬ 
cate  that  there  must  be  a  source  at  depth  in  the  upper 
elevations  of  Caribou  Peak.  Yet  the  borehole  location  is 
only  60  m  in  elevation  below  the  Peak’s  crest.  The  fact 
that  multiple  fractures  exist,  each  exhibiting  reflectors 
corresponding  to  water  flow,  indicates  that  the  source  is 
probably  not  a  single  spring  outcropping  near  the  crest, 
but  rather  a  more  pervasive  aquifer.  It  should  be  noted 
that  no  springs  have  been  observed  within  this  area.  The 
ground  temperatures  also  require  that  water  moving 
through  the  fractures  be  sufficiently  warm  to  maintain 
a  localized  unfrozen  condition  within  the  otherwise 
permanently  frozen  bedrock  that  encompasses  them.  In 
addition,  the  chemical  and  NMR  analyses  indicate  the 
source  must  be  fresh  water. 

Our  hypothesis  concerning  the  origin  of  the  ground- 
water  is  based  upon  these  observations  and  the  geologic 
limitations  of  this  site.  The  general  north-south  strike 
and  low  angle,  easterly  dip  of  the  fractures  indicate  that 
they  bisect  the  hill  and  probably  continue  through 
Caribou  Peak  to  its  southern  and  southeastern  slopes 
(Fig.  14).  Vegetative  patterns  and  associated  rock  out¬ 
crops  identified  on  aerial  photographs  also  appear 
continuous  from  the  north  to  south  face.  If  they  are 
continuous  through  the  hill,  the  fractures  must  intersect 
the  non-permafrost  zone  below  the  south  face.  An 
aquifer  within  the  unfrozen  bedrock  here  could  provide 
a  source  of  water  for  the  fractures. 

Borehole  temperature  measurements  at  mid-slope  on 
the  southern  face  confirm  that  permafrost  is  absent.  Mean 
annual  surface  temperatures  range  from  0.4  to  2.2°C,  and 
ground  temperatures  at  10  m  depth  range  from  0.9  to 
1.9°C  (Collins  et  al.  1988).  Since  groundwater  tempera¬ 
tures  approximate  the  mean  annual  air  temperatures,  we 
anticipate  that  groundwater  within  this  non-permafrost 
area  would  have  a  temperature  somewhere  within  the 
range  of  0.9  to  1 .9°C.  This  range  would  be  consistent  with 
groundwater  temperatures  in  the  Fairbanks  region  of  an 
average  1.7°C  as  reported  by  Williams  (1970). 
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Figure  13.  Idealized  diagram  showing  a  possible  relationship  between  water¬ 
bearing  fractures  and  vegetative  lineations. 
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Thus,  if  we  assume  that  the  fractures  intersect  an 
aquifer,  water  may  be  transmitted  within  them  and 
along  their  entire  length  from  the  warmer  southern  side 
to  the  northern  face  (Fig.  14).  Discharge  along  the 
length  of  the  fracture’s  intersection  with  the  colluvium 
would  be  consistent  with  the  continuous  nature  of  the 
vegetative  lineations.  Water  can  move  within  the  frac¬ 
tures,  however,  only  if  the  combination  of  permafrost 
thickness,  permafrost  temperatures,  groundwater  tem¬ 
peratures  and  flow  zone  properties  result  in  sufficient 
flow  rates  and  heat  flux  to  maintain  an  unfrozen  zone 
continuously.  The  water’s  temperature  must  remain 
above  freezing  while  it  is  flowing  through  the  frozen 
rock. 

The  air  and  ground  temperature  measurements 
suggest  that  the  perennially  frozen  bedrock  probably 
does  not  exceed  -1.3°C  in  temperature  at  its  coldest 
point  near  the  10-m  depth.  The  thickness  of  permafrost 
on  the  north  face  is  unknown,  but  it  probably  is  not  very 
thick,  as  is  characteristic  of  this  region  (Ferrians  et  al. 
1969,  Williams  1970).  An  approximate  depth  to  the 
base  of  permanently  frozen  ground  can  be  determined 
from  the  ground  temperatures  in  boreholes  K20A  and 
K2,  if  we  assume  that  q  and  K  (eq  4  and  5)  are  constant 
at  this  location.  Calculated  mean  annual  temperature 
gradients  yield  an  estimated  depth  to  the  base  of  per¬ 
mafrost  (0°C  depth)  of  1 0  m  at  K2  and  40  m  at  K20A. 
These  relatively  shallow  depths  indicate  a  thinning  of 
perennially  frozen  ground  to  the  crest  of  Caribou  Peak, 
as  expected. 

If  groundwater  is  moving  within  the  fractures  from 
the  southern  to  northern  side,  we  can  assume  that  its 
initial  temperature  would  be  approximately  1 .5°C  and 
that  it  must  flow  through  the  estimated  35  m  of  frozen 
materials  in  the  K20  borehole  area,  cooling  to  slightly 
above  0°C.  The  frozen  ground  temperatures  would 
gradually  decrease  from  O^C  at  about  40  m  depth  below 
the  surface  to  about- 1 .3°C  at  1 0  m  depth  beneath  it.  The 
unfrozen  zone  or  talik  adjacent  to  the  fracture  would 
extend  outward  and  parallel  to  the  fracture  a  short 
distance,  dependent  upon  the  groundwater  heat  flux 
(Fig.  14). 

Once  a  flow  system  is  established,  a  thermal  equilib¬ 
rium  would  necessarily  develop,  offset  only  by  seasonal 
warmings  from  groundwater  recharge  by  major  pre¬ 
cipitation  events  during  summer.  Such  seasonal 
warmings  were  suggested  by  the  slight  increases  in 
borehole  K20 ground  temperatures  from  May  to  August. 
Heat  conducted  by  the  water  from  the  southern,  non¬ 
permafrost  zone  must  maintain  the  unfrozen  condition, 
or  the  zone  would  refreeze,  with  ice  formation  terminat¬ 
ing  groundwater  movement.  Once  such  a  system  is 
established  and  stabilized,  heat  conducted  by  the 
groundwater  could  maintain  the  unfrozen  state  (e.g.. 


Williams  1970).  If,  however,  groundwater  discharge 
through  the  zone  decreases,  temperatures  at  the  outer 
fringe  of  the  zone  would  begin  to  decrease  and  initiate 
freezing-in  of  the  fracture  talik. 

Rigorous  calculations  to  determine  if  this  is  theoreti¬ 
cally  possible  at  this  site  would  require  data  not  avail¬ 
able  to  us,  principally  the  actual  location  of  the  aquifer 
relative  to  fracture  discharge,  the  relative  “permeabil¬ 
ity”  of  the  fracture  that  controls  the  rate  of  flow,  and  the 
hydraulic  head.  Movement  within  fractures  can  be 
relatively  rapid,  far  exceeding  that  which  would  depend 
upon  the  permeability  of  the  bedrock  itself.  Groundwa¬ 
ter  movement  from  subpermafrost  aquifers  through  a 
permafrost  horizon  in  bedrock  fractures  has  been  de¬ 
scribed  at  various  locations  by  Hopkins  et  al.  (1955), 
Williams  and  Van  Everdingen  (1973),  Tolstikhin  and 
Tolstikhin  (1974)  and  Ivanik  and  Veselovskiy  (1976), 
among  others.  Physical  reasoning  suggests  that  flow 
could  occur  at  a  rate  sufficient  to  maintain  these  near 
0°C  water  temperatures,  inhibit  freezing  and  maintain 
the  locally  thawed  zone  through  this  relatively  thin, 
warm  permafrost  in  a  manner  analogous  to  water  move¬ 
ment  through  fractures  in  temperate  glacier  ice. 


CONCLUSION 

The  radar  survey  has  identified  a  fracture  in  the 
quartz-muscovite  schist  bedrock  with  groundwater 
flowing  in  it  as  the  source  of  the  anomalous  unfrozen 
zone  within  borehole  K20.  The  surveys  indicate  that 
several  fractures  with  groundwater  flow  exist  in  the 
surveyed  area  of  the  slope,  but  were  not  intersected  by 
other  drill  holes.  The  uniformity  of  the  near  0°C  ground 
temperatures  measured  over  the  last  3  years  indicates 
both  the  permanence  of  the  unfrozen  zone  and  a  require¬ 
ment  for  the  flow  of  water  through  the  fracture.  The 
distribution  of  the  water-bearing  fractures  in  relation  to 
the  vegetative  pattern  of  the  north-facing  slope  suggests 
that  they  act  as  both  water  sources  and  water  transmitters 
from  the  peak  to  the  valley  bottom,  enhancing  the 
growth  of  vegetation  where  they  outcrop  below  the  col¬ 
luvial  cover.  The  ultimate  source  of  the  groundwater, 
which  constantly  remains  above  freezing  and  thus  re¬ 
mains  an  active  flow  system  annually,  is  unknown,  but 
may  include  an  aquifer  within  the  hill  that  penetrates  the 
permafrost  beneath  the  northern  slope  along  permeable, 
hydraulic  anisotropies  created  by  the  bedrock  fractures 
that  intersect  the  unfrozen  south-facing  slope  (Fig.  14). 
Future  work  to  verify  the  source  of  water  in  the  fractures 
should  include  detailed  field  studies  of  the  groundwater 
system  through  radar  surveying  of  the  fracture  distri¬ 
bution  in  Caribou  Peak,  and  drilling  and  analysis  of  the 
groundwater  flow  pattern  and  temperatures. 
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a.  Cross  section  of  Caribou  Peak,  view  to  east 
( top)  and  map  view  of  horizontal  slice  through 
Caribou  Peak. 
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b.  Cross  section  parallel  to  transect  I ,  view  to  north. 

Figure  14.  Schematic  model  of  the  apparent  relationship  beween 
fractures,  groundwater  occurrence  and  vegetation.  The  hypothesized 
source  of  groundwater  is  an  aquifer  beneath  the  non-permafrost  southern 
slope  that  is  intersected  by  the  fractures,  which  then  transmit  this  water 
through  the  permafrost  horizon  of  the  northern  face. 


The  presence  of  unfrozen  zones  of  this  type  requires 
consideration  when  predicting  the  distribution  of  per¬ 
mafrost  in  the  Interior  region.  Similar  patterns  to  the 
vegetation  on  other  north-facing  slopes  suggest  that 
thawed  zones  and  groundwater  systems  associated  with 
bedrock  fractures  may  occur  elsewhere  within  the  wa¬ 
tershed  region. 

Groundwater  flow  on  this  northern  slope  discharges 
into  Poker  Creek.  This  discharge  contributes  to  the 
baseflow  of  the  watershed,  but  has  not  previously  been 
considered.  It  may  be  important  to  understanding  higher 
than  normal  winter  flows  and  baseflow  recession  in 
watersheds  underlain  by  discontinuous  permafrost  (e.g.. 


Dingman  1971.  Haugen  et  al.  1982),  as  well  as  the 
source  of  water  forming  extensive  aufeis  within  the 
Poker  Creek  valley  bottom  each  winter. 

Finally,  this  study  has  verified  the  ability  of  impulse 
radar  operating  at  high  frequency  to  identify  localized 
unfrozen  zones  associated  with  groundwater  flow  in 
perennially  frozen  bedrock. 
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in  the  direction  of  the  slope,  trending  generally  N20°E  and  dipping  about  9°  east.  Geologic  logs  indicate  that  the  drill  hole  intersected  a  fracture 
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that  the  water  was  fresh  and  not  highly  mineralized.  Vegetation  patterns,  coupled  with  the  borehole  location  and  fracture  orientation,  suggest 
that  flow  originates  within  the  upper  and  central  parts  of  the  peak  and  discharges  into  the  valley  of  Poker  Creek.  The  source  of  the  groundwater 
is  unknown,  but  appears  to  be  an  aquifer  in  the  south-facing,  non-permafrost  side  of  Caribou  Peak  that  is  intersected  by  the  north-south  striking 
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